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Abstract: Responding mechanically to environmental requests, muscles show a surpris- 
ingly large variety of functions. The studies of in vivo cycling muscles qualified skeletal 
muscles into four principal locomotor patterns: motor, brake, strut, and spring. While 
much effort of has been done in searching for muscle design patterns, no fundamental con- 
cepts underlying empirically established patterns were revealed. In this interdisciplinary 
study, continuum mechanics is applied to the problem of muscle structure in relation to 
function. The ability of a powering muscle, treated as a homogenous solid organ, tuned to 
efficient locomotion via the natural frequency is illuminated through the non-linear elastic 
muscle moduli controlled by contraction velocity. The exploration of the elastic force pat- 
terns known in solid state physics incorporated in activated skeletal and cardiac muscles 
via the mechanical similarity principle yields analytical rationalization for locomotor muscle 
patterns. Besides the explanation of the origin of muscle allometric exponents observed 
for muscles in legs of running animals and wings of flying birds, the striated muscles are 
patterned through primary and secondary activities expected to be useful in designing of 
artificial muscles and modeling living and extinct animals. 

PACS: 89.75.Kd, 89.75.Da, 87.10.Pq, 87.19.Ff 
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I. INTRODUCTION 



The mechanical role of muscles varies widely with their architecture and activation con- 
ditions. Striated (skeletal and cardiac) muscles are diverse in their contractive interspecific 
and intraspecific functional properties observed among and within animal species, neverthe- 
less, "the smaller muscles and muscles of smaller animals are quicker". This generic feature 
of skeletal muscles was established by Hill [1]. More recently, the physiological adaptation 
of skeletal muscles resulting in beneficial changes in muscle function has been recognized by 
a number of investigators. It was learned that long-fibre muscles commonly contract at over 
larger length ranges and relatively higher velocities producing the greatest muscle forces 
the lowest relative energetic costs [2]. Muscles having shorter fibres expose smaller length 
change, but their cost of force generation is relatively less, e.g. [3]. Searching for determi- 
nants of evolution of shape, size, and force output of cardiac and skeletal muscle, a little 
is known about the regulation of directional processes of mass distribution [4,5]. Although 
skeletal muscles grow in length as the bones grow, most studies only involve force increasing 
with cross-sectional area. Following the idea that the muscle force production function is a 
critical evolutionary determinant [5], I develop a physical study of muscle form adaptation 
to a certain primary activity with growth of size preserving muscle shape. 

When designing architecture of the striated muscle built from repeating units (fibres 
and sarcomeres) at least three distinct muscle activities should be distinguished [5]: (i) 
the concentric contraction defined as the production of active tension while the muscle 
is shortening and performing positive work, (ii) the eccentric contraction defined as the 
contraction during lengthening performing negative work in a controlled fashion, and (iii) 
the isometric contraction when the muscle force output is produced without changing of 
length and performing net work. The corresponding mechanical work patterns called by 
Russel et al. [5] as "concentric work" and "eccentric work" (that might be extended by 
"isometric work") were carefully studied via in vivo measurements of length-force cycling of 
individual skeletal muscles in active animals, such as (i) the pectoralis in flying birds, (ii) 
leg extensors in running cockroaches, and (iii) gastrocnemius in the level running turkey. 
The corresponding muscle locomotor patterns were called as (i) motor, (ii) brake, and (iii) 
strut functions [6]. 

The seminal research by Hill [1] on dynamics of electrically stimulated isolated muscles 
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was restricted to a single isotonic shortening. The studies of the corresponding motor func- 
tion resulted in famous force-inverse-velocity master curve remaining the major dynamic 
constraint of all real (slow-fibre, fast-fibre, and superfast) muscles [7] and computationally 
modeled muscles, e.g. [8]. Besides, other two fundamental rules of muscle dynamics were 
noted by Hill [1]. Examining hovering humming and sparrow birds, he recognized that the 
"frequencies of wings are roughly in inverse proportion to the cube roots of the weights, 
i. e. to linear size" . Moreover, because the linear proportionality between the stroke period 
and body length was found equally in electrically stimulated isolated muscles, the intrinsic 
frequency- length feature constrained by scaling rule T m l oc L" 1 beyond the nervous control 
is likely more universal than previously appreciated. Second velocity-length Hill's constraint 
states that "the intrinsic speed of muscle has to vary inversely to length", i.e. V m oc L" 1 . 
Both Hill's scaling rules remain a challenge to viscoelastic models of transient-state mechan- 
ics and other existing theories of muscle contraction [9]. 

The earliest theories of muscle motor function supposed muscle to be an elastic body 
which, when stimulated, was converted in an active state containing elastic energy causing 
the muscle to shorten. Such elastic- energy theories failed to explain mechanisms of the force 
production in terms of viscoelastic characteristics. To a certain extent, poor experimental 
approaches providing often conflicting clues to muscle structure in relation to function may 
explain a little progress in understanding of contractile properties of a muscle [9,4] . Moreover, 
physiological muscle properties accounted for theories of muscle contractions developed at 
both molecular and macroscopic scales are primarily focused on the reproduction of force- 
velocity curve [9]. Besides, the existing phenomenological frameworks such as Hill-type 
muscle models only mimic the proper mechanical characteristics of muscles by means of 
passive viscoelastic springs attached to muscle contractive element in series [10,11,3] or in 
parallel [12] and recruited when muscle is activated. By ignoring the proper muscle function 
of force production and force transmission throughout the muscle organ, these models are 
able to explain no one of Hill's principal constraints in muscle dynamics. On the other hand, 
there exist experimental evidences of the adaptive ability of skeletal muscle to exchange 
elastic strain energy during force production [2]. In line with this concept, it has recently 
communicated on a possibility of the modeling of the adaptive muscle elasticity by elastic 
force patterns [13]. 

In the present paper, I develop an integrative theoretical framework to the problem of 



3 



forces, structure, and contractive non-linear dynamics in striated muscles. Instead of Hill- 
type modeling of in vitro motor function, e.g. [3], brake function, e.g. [12, 2], and strut 
function, e.g. [14], or study of muscle design by means of simulation of phenomenological 
force-length and/or force- velocity constraints [8], the powerful method of continuum me- 
chanics generally providing macroscopic characterization and modeling of soft tissues, e.g. 
[15, 16], is employed. By further exploration of the elastic force patterns, I propose a self- 
consistent depiction of the three velocity-distinct characteristic points well distinguished in 
all in vivo force-length loops of the naturally activated muscles. Unlike the earliest elastic 
theories based on minimization of energy, I develop the physical concept of similarity between 
the force output and reaction active elastic forces that permits to avoid the details of muscle 
activation process. The theory is validated by a comparison to phenomenological scaling 
rules including both mentioned Hill's dynamic constraints and therefore may be hopefully 
helpful in designing artificial muscles [15] and modeling living and extinct organisms [17]. 

II. THEORY 

A. Theoretical Background 

1. McMahon's scaling models 

The engineering models by McMahon [18, 19] develop previous Hill's approach to the 
problem of scaling of parameters of animal performance to body weight W = Mg. Using Hill's 
geometric similarity models [1,19] equally applied to animal body, long bone, or individual 
muscle, each one was approximated by a cylinder of longitudinal length L and cross-sectional 
area A (or diameter D ^ \/A). Then, the assumption on the weight-invariance of for the 
tissue density, namely 



was adopted. In mammalian long-bone allometry, this invariant was verified and observed 
with a high precision [20]. Mechanical models of bending bones and shortening muscles 
were introduced by McMahon via the weight-invariant elastic modulus E tiss , stress o- tiss , 
and strain e tiss , namely 
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Here AL (= L — Lq) is the length change accompanied by the force change AF (= F — Fq) 
counted off from the resting length Lq. 

While searching for functional mechanical patterns of biological systems determined by 
maximal forces using Eqs. (CQ) and (|2J), the maximal stress/strain scaling relations 

« ™ X, \ oc WV\ and a W/ 1 / 5 , (3) 

could be readily derived from McMahon's geometric (isometric), elastic (buckling stress) and 
static (bending elastic stress) similarity models distinguished through McMahon's scaling 
relations 

L g eom ^ D, L elast cc D 2/3 , and L stat cc D 1/2 . (4) 
Instead, the maximum stress and strain 

'S? oc a W°, (5) 

were postulated (see Table 4 in [19]) extending groundlessness his exact result for the mean 
stress c^Iast oc W°, obtained within the static stress similarity model (see Fig. 1 in [19]). 
The improved self-consistent maximal stresses shown in Eqs. (|2D follow straightforwardly 
from McMahon's cross-sectional 

AgT> a W^, a W 3/4 , and A<£g a H/ 4 / 5 (6) 

applied to Eq. (j2J), along with McMahon's idea on the dominating gravitational forces in 
bones, muscles, and bodies, i.e. AF ^ gM b ^ gM m ^ VK. As shown in [20], the structure 
of long bones is driven by peak muscle forces, but not by gravity. 



2. Muscle shape and structure 

After Alexander [21], the physiologic cross-sectional area A 0m (PCSA) of the isolated 
skeletal muscle m of mass M m composed of N bundles of masses was commonly esti- 
mated, e.g. [22], with the help of the cylinder-geometry relation Ai = rrii/ p musc Li, where 
Li is directly measured muscle fibre length. The spindle-like shape of the muscle as whole 
organ was therefore determined by the muscle PCSA, namely 

-sr-^ M m 1 1 v-^ m « \ 

A 0m = 2_^ A i = j — ' and 7 — = 7p2^T"' w m = Z^ mi > (') 

1=1 i=l i=l 
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resulted in the sum of areas of muscle and the muscle length Lo m of the parallel-linked 
contractible subunits described statistically by the length-unversed sum weighed by masses. 
Such a coarse-grained characterization of the muscle structure generally ignores the arrange- 
ment of muscle fibres relative to generated force axis, distinguished by pinnate angles. 

In scaling models, the evolution of the muscle structures across different-sized animals of 
body mass M is observed statistically via allometric exponents a m , l m , and a m determined 
by common rules [21,23,25]: 

A 0m a M£» , L 0m a M l ™ , and M m cx M 1+a ™ , (8) 

where the muscle mass index a m plays the same role as Prangel's index j3 in bones, as noted 
in [26]. 

When the muscle-density invariance employed implicitly in Eq. (j7j) and specified in 
Eq. (TjQ) is applied to different skeletal muscles, the muscle shape approximated by cylinder 
geometry is also preserved. Consequently, the muscle functional volume 

A m L m = A 0m L 0m = — -, with p 0m = p musc cx M° t a M°, (9) 

POm 

holding in all muscle work loops plays the role of the mechanical muscle invariant. This 
condition is ensured by functional change Ap musc / p musc not exceeding 5% [24]. Hence, the 
function-independent muscle-shape constraint [13] 

a m + l m = l + a m (10) 

straightforwardly follows from Eqs. (JSJ) and Q. Likewise the case of hindlimb mammalian 
bones of the mean structure a[ exp) = 2dJ; exp) = 0.752, l ( ° xp) = 0.298, and (3^ = 0.04 [20,26], 
Eq. (TlDl is also empirically observable in muscle allometry (see analysis in Table 5 below). 

B. General Muscle Characterization 

1 . Maximal force and stress 

In in vivo work loops, the muscle locomotor patterns can be generally specified regardless 
of details of activation-deactivation conditions. In Fig. 1, the linear-slope characteristics 
L\ m can be introduced in the force-length cycling by the domains: L 2m < L\ m < L 3m w L 0m , 
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for the motor function, L 2m > L\ m > L 3m , for the brake function, and by L 2m > L\ m > 
pa Lo m , for the strut function showing nearly isometric muscle contractions. 

Place Fig. 1 

Moreover, such a qualitative general characterization of the activated individual muscle 
m of resting length Lo m can be rationalized on the basis of common two-point force-length 
characterization, namely 

(Lim) = Fj%$ = F 2m and F^l (L lm ) = F Xm , (11) 

introduced by the maximum force F 2 ™ and the optimum muscle length [28, 29] L\ m . 
The instant dynamic length L m = L lm ± AL lm is counted off from the characteristic point 
Li m via the optimum length change ALi m (I lip shown in Fig. 1 for all functions. 

First, the linearization of the in vivo muscle force- length curve allows one to determine 
trial peak stress and corresponding strain by 

p(max) a j (max) 
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The corresponding force change AF^l^ observed near the optimum force F^^ (ITT]) pro- 
vides 

p(max) _ p(cxp), T \ , a u(max) _ p , Tz(mia) a r (max) /19\ 

1 muse ~ 1 muse K^lm) ~ L -^ 1 muse lm i ly musc L -^ 1J lm \ LO ) 



that determinates effective muscle stiffness and effective modulus, respectively 
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and iff - ^ = (14) 
following from Eqs. ffT2l) and ([TBI) . 



2. Active stiffness and resonant muscle mechanics 

Secondly, treating the maximum-force crossover state (fTT|) as the generic transient-neutral 
state [26] the resonant frequency 1/T^nu^ = T 2 ~ m ~ related to point 2 in Fig. 1 associated with 
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maximum efficiency of muscle cycling, e.g. [29], can also be introduced as natural frequency 
[19], namely 

pr? / p(max) / « p (max) \ X / 2 . 

m— 1 M v 2m -L^musc I ^-T muse \ J- /1 _n 

r 2m ~27rW— -~y— r(max) 7—, (15) 

V M m V pom \ Fi usc ; / ^2m 



and analyzed by Eqs. (j9]) and (p2 

One can see that Eq. (fl5|) yields first Hill's general constraint discussed in Introduction. 
However, the following three conditions are required: (i) the preservation of dynamic func- 
tional volume (j2J), (ii) the weight-invariance of the elastic modulus Emusc (ED, and (iii) the 
existence of force similarity between the exerted force F^nu^c and its change AFmusc ( |T3j) . 
Therefore, the muscle force- similarity principle, namely 

Fmusc = AF musc = Fp ro d = F e i as i = AF e i a sti (16) 

implying a coexistence of all forces in biomechanically equivalent states [26] must be adopted. 
Here the active elastic force AF e i ast (shown schematically as F act in Fig. ID) is also included. 
The total state-transient elastic force F e i ast is the superposition of common passive elastic 
force F pass provoked by external loads and active elastic force AF e i ast caused by the pro- 
duction force Fp ro d. The correspondence sign = indicates that though the involved physical 
characteristics belong to the same mechanical state, they may differ in both physical and 
numerical parameters stipulating this state. 

Given that the peak active muscle stress a m always exceeds the corresponding passive 
stress, e.g. [14], in further I focus on the fully activated transient states described by 

AF e i as t ,-, AL m 

cr m = = E m — — . (17) 

Unlike Eqs. ([2]) and (1121) . o~ m is the true intrinsic elastic stress in a certain (not specified) 
dynamic state. This reveals the maximum-amplitude elastic force of the fully activated 
muscle 

AF e i as t = AF m = K m AL m = E m A m m (18) 



and in turn provides the corresponding active muscle stiffness 

K m = E m ^. (19) 

The underlying mechanical sarcomere elastic stiffness K s is related via the muscle-volume 
average, namely 

K m = — ^— [ K s {r m ) d 3 r m , (20) 



originated from end-to-end intercellular overlapping [31, 12]. 
The muscle energy change 
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stored or released during active-period contraction provides the mechanical cost of energy 
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These relations demonstrate how the observable mechanical characteristics can be linked to 
the underlying muscle elastic forces using the force-similarity principle (fT6l) . In turn, the 
contraction velocity 
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is defined by the instant velocity V m (t) averaged over activation time At m . 



3. Fast and slow activated muscles 



According to the most general classification of diverse muscles, three types are conven- 
tionally distinguished: red (slow fibre) muscles, white (fast fibre) muscles, and intermediate 
type, mixed fibre muscles. Although collective mechanisms of muscle contractions are poor 
understood, e.g. [32], physically, the two limiting situations of dynamic accommodation 
of local forces generated by cross bridge attachments can be generally rationalized. As 
schematically drawn in Fig. ID, in an activated muscle, the dynamic process of equilibra- 
tion between the production intrinsic forces and external loads (not shown) is followed by 
the spatiotemporal relaxation of elastic forces. For the simplest case of slow muscles, the 
dynamic equilibration occurs via the slow channel of relaxation, assumably common for 
both active, F^ r l °^\ and passive elastic forces. Since passive forces in solids are short of 
range [33], both the forces are proportional to muscle surface. In contrast, it is plausible 
to adopt that in fast muscles the fast-twitch fibres transmit the locally generated forces in 
all directions, i.e. along and across fibres, resulting in the overall maximum force output 
Fplod^ t° be linear with dynamic muscle volume. Basing on such a general physical picture, 
a function-independent and regime- independent characterization of the force production 
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function, namely 



Fplof a ArmLrm and cc A rm , with r = 1,2, and 3, (24) 

is proposed via the force-size scaling rules for all three distinct states shown in Fig. 1 and 
hereafter distinguished by symbol r. 

The widely adopted by biologists linear-displacement regime is discussed in Eq. (j2J) 
via AL cc L resulted in the weight-independent strain (jSJ). The corresponding optimum- 
velocity regime r = 1, attributed to the instant length-independent elastic strains, etf'' = 
\L m — L 3m \/L m cc L Q m (jSJ) with L m lying between L lm and L 3m rs L 0m , is now clarified by 
the scaling equations 

characteristic of fast and slow muscles. Such a muscle description follows from the similarity 
(flE]) between the active elastic force AF lm = AF^J = Ei m Ai m Ei m (fT8l) and corresponding 
production force f l24l) . The optimum force- velocity muscle mechanics is rationalized below 
in Table 1 and then tested by empirical data. 

Similarly, the bilinear-displacement regime r = 2 introduced by the dynamic length 
change AL 2m = \L m — L\ m \ cc L 2 m , with L m lying between L 2m and L\ m , and the maximum 
active elastic force AF 2m = AF^f 1 = E^m^m^m^ results in the maximal elastic 
moduli 

p(max) _ Tp(fast) r j p(max) _ rp{slow) „ r -1 / 0K \ 

adjusted with the muscle production function f[2~4"j) via the force similarity principle (TTB1) . 
Finally, the high-velocity trilinear regime r = 3 is suggested by the moderate-force and 
moderate-elastic muscle determined by 

E ( p d) = cc L- 1 and E^f = cc L~ 2 . (27) 

This condition specifies point 3 in Fig. 1, along with the underlying cubic-power muscle 
displacements AL 3m cc scaled by dynamic L m lying above or below the characteristic 
length L 3m in any muscle acting as motor, brake or strut (see Fig. 1). 
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C. Muscle Functions 

Likewise the naturally curved mammalian long bones biomechanically adapted to the 
maximum longitudinally bending [20, 26], the muscle motor function is assigned to locomotor 
muscles showing concentric positive work exerted by elastic bending forces. Given that the 
elastic force patterns coincide for bending and torsion [26], both kinds of unpinnate and 
uni-pinnate skeletal muscles, having respectively close to zero and non-zero fixed pinnate 
angles, may be expected to be structured by the same motor function. The specific-function 
mechanical characterization is described in Appendix B and results are summarized in Table 
2. 

III. RESULTS 

A. Assumptions and predictions 

The following assumptions are made regarding elastic striated muscles: 

1. The muscles are considered at macroscopic scale as individual homogeneous organs. 
Within the continuum mechanics, the coarse-grained approach ignores the details of hetero- 
geneous microstructure and pinnate angles. 

2. When activated under different boundary loaded conditions, the muscles do not un- 
dergo changes in shape and whole volume. The emerging elastic fields follow patterns es- 
tablished for long solid cylinders. 

3. The mechanical similarity adopted between the extrinsic forces exerted by the muscle 
and intrinsic elastic reaction forces, as well as the dynamic similarity adopted for contraction 
velocities and frequencies are observable in all biomechanically equivalent states. 

4. The natural ability of the non-linear elastic tuning of fast and slow muscles to distinct 
locomotor states can be characterized by the elastic moduli sensitive to evolving dynamic 
variable associated here with the regime-characteristic muscle length. 

The function-independent mechanical characterization of muscles is provided in Table 1. 

Place Table 1. 
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The specific case of muscle structure accommodation in the bilinear regime is described 
in Table 2. 

Place Table 2 

The rules of mass distribution across and along the muscle axis provided in Table 2 in 
terms of the muscle-structure scaling exponents [dam, hm] are characteristic for slow, fast, 
and mixed muscles producing maximum force. In Table 3, these scaling rules are compared 
with the finding for the optimal- force state [ai m , li m ] and moderate- force state [a^m, hm}. 

Place Table 3 

The dynamic characteristics of distinct-velocity contractions are predicted in Table 4. 
Place Table 4 

The consequences of the theoretical scaling framework are: 

1. The peak forces generated in all regimes scales as muscle volume or surface in fast or 
slow muscles, respectively. 

2. A general, function-independent mechanical description of the striated muscle acti- 
vated in the liner-displacement regime is predicted for each type of muscles (Table 1). 

3. The muscle- type independent locomotor functions and related mechanical and dynamic 
characteristics of the striated muscle activated in the bilinear regime are predicted (Table 
2). 

4. The muscle-type independent varied dynamic structures are predicted for all muscle 
regimes and functions (Table 3). 

5. The function- independent dynamic scaling characteristics are obtained in Table 4 for 
all type of muscles. 

In what follows, all theoretical findings are tested by the available from the literature 
data. 
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IV. DISCUSSION 



"What determines the shape, size, and force output of cardiac and skeletal muscle?" 
(Louis Sullivan quoted in [5]). The provided coarse-grained study of conservative stri- 
ated muscles suggests that the size- dependent peak elastic forces determine fiber-type- 
independent patterns of the functionally adapted structures preserving muscle shape. The 
size-dependent peak force output is determined by the muscle volume and area for white 
and red muscles, regardless of muscle structure and function. 

A. Function against structure 

1. General muscle characterization 

Being composed of bundles of muscle fibres including all other contractible components 
(neural, vascular, and collagenous reticulum), the striated muscle is thought of as a heteroge- 
neous continuum medium transmitting the produced tension internally and externally, e.g. 
[34]. Primarily, I address the problem of mechanical design of striated muscle to a general, 
function-independent characterization of the individual muscle organ loaded by tension, re- 
action, and gravity through tendons, ligaments, and bones. My non-energetic approach is 
physically grounded by the existence of linear force- displacement regions (shown by the solid 
arrows in Figs. 1A, IB, and 1C) in all in vivo work loops regardless of dynamic details of 
approaching to the maximum exerted force F^!^ . Hence, the mechanical characterization 
of the maximum-force activated muscle arises from the muscle stiffness (1141) underlaid 

by sarcomere stiffness A"i max ' ) (1201) . Consequently, all forces involved in muscle contraction 
following by active and passive elastic strains allow common mechanical description (shown 
in Fig. ID) not depending on their biochemical, inertial, or reaction origin. 

The analytical justification of Hill's first frequency-length constraint arises from the analy- 
sis of Eq. (1151) that requires eventually the usage of the similarity between all intrinsic muscle 
forces, Eq. (|T6|) . The constraint T m x oc L" 1 and other mechanical characteristics for slow 
muscles accumulated in Table can be applied to steady- speed regimes of locomotion modes 
where all forces are generally equilibrated and controlled by slow-fibre muscles [35]. In the 
case of non-steady transient locomotion when fast-twitch fibres and nervous control are ad- 
ditionally requested [35], Hill's first constraint transforms [by Eqs. (|T5l) and (1251) ] into a new 
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one, X' 1 oc L m ' oc 1 /V^ t (Tables 1 and 4), well known for animals running with maximal 
optimum speed [36,37] Ktm^ ^ ^/aft oc a V-^m- We have therefore demonstrated how 
the suggested dynamic similarity establishes a link between the body-propulsion speed and 
locomotor-muscle contraction velocity, also described by Rome et al. [38]. Being united with 
the muscle-force similarity, both constraints yield mechanical similarity, the key principle 
explored in this research. 

2. Maximum force output against structure and velocity 

In muscle physiology, the functional effect of muscle conceptual architecture simply states 
that muscle force output is proportional to PCSA. The proposed study of adaptation of the 
muscle structure via the force production function seems to be in qualitative agreement with 
this statement, because in all cases exposed in Eq. ( 1241) the muscle force output is propor- 
tional to A m . Such a simplified treatment of the fast-muscle mechanics (formally substituted 
by that for slow muscles) arrived at the widely adopted opinion that the peak muscle stress 
Fprod ' ] /Am > specifying the case of slow muscles in linear dynamic regimes with ffm™' a L? m 
(Table 1), is generic for any muscle, as already discussed in Eq. (jSJ). Although the proposal 
on scaling of the maximum production force (and active stress) with muscle size (124"]) is a 
challenge for further research, the provided fairly general physical grounds are supported 
by empirical observations by Marden and Allen [39]. They established statistically that the 
maximum force output in all biological (and human-made) motors falls into two fundamental 
scaling laws: (i) in fast-cycling motors, presented by flying insects, bats and birds, swim- 
ming fishes, and running animals it scales as [motor mass) 1 and (ii) in slow-cycling motors, 
such as myosin molecules, muscle cells, and some (unspecified) "whole muscles" the force at 
output scales as (motor mass) 2 / 3 . The "motor mass" was associated with muscle (and fuel) 
mass. That fact that the authors observed muscle motors from sarcomere to whole muscle 
organ passing through the single-fibre level of muscle organization, makes a basis for the 
discussed below micro-macro scale correspondence. 

The proposed treatment of the in vivo force-length curves is provided for three distinct 
force- velocity characteristic points (shown in Fig. 1) correlated by the inequalities 

F2m > F\ m > F 3m and V 2m < V\ m < V 3m . (28) 
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These three generic function-independent states are associated with the linear (r = 1), 
bilinear (r = 2), and trilinear (r = 3) muscle dynamics determined via the muscle elastic 
moduli E rm in Eqs. (1251) . (12B")) . and (127)) . respectively. The mechanical characterization 
of slow and fast striated muscles is therefore provided in terms of the maximum (AF 2m ), 
optimum (AF lm ) and moderate (AF 3m ) active elastic forces developed at the measurable 
maximum (V^m), optimum (Vi m ), and moderate (V2 m ) contraction velocities (Table 4). The 
stabilization of the dynamic regimes is expected at the natural frequencies, which also are 
scaled in Table 4 to the dynamic length L rm . 

3. Muscle functions against size and shape 

Searching for answer on "what features make a muscular system well-adapted to a specific 
function?" [28], it has been shown preliminary [13] that such features are related to natural 
conditions of the stabilization or tuning to the moderate- velocity regime r = 2 via the mean 
dynamic length of the fast-twitch fibers adapted by the best way to one of the patterns of 
muscle locomotor functions. In this study such features specify the role of slow-twitch fibers. 

The elastic-force patterns underlying concentric, eccentric, isometric, and cardiac con- 
tractions are suggested in Eqs. (1351) . (1391) . (1421) . and (144)) . respectively. The solutions to the 
muscle-force and muscle-shape constraints are accumulated in Table 2 as patterned functions 
well distinguished by the muscle structure parameter (7] m = din A m /d In L m ) established for 
the motor (rji = 4), brake (rj 2 = 3), strut (773 = 00) skeletal muscles, an extended by the 
pump (rj 5 = 1) cardiac muscle and one spring = 2) striated muscle. These structurally 
adapted muscles are thought of as to be suited to efficient work during powering when, 
respectively, shortening (m = 1), lengthening (m = 2), or remaining in the nearly isometric 
dynamic state (m = 3), high-pressure-resistant state (m = 4), and likely energy-saving state 
(m = 5). The found new pump function is in accord with the observation by Russel et 
al. [5] that "the heart chamber, unlike skeletal muscles, can extend in both longitudinal 
and transverse directions, and cardiac cells can grow in length and width", that implies 
V5 < Vii^i or rj 3 . Given that only a few patterns exist in elastic theory of solids [26], it is 
not striking that the spring, brake, and motor functions resembles McMahon's "geometric", 
"elastic", and "static" stress similarities discussed in Eqs. ([3]) and (jl)). 

In Table 3, conceivable stable dynamic structures corresponding to muscle activity in 
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different dynamic regimes are analyzed. As in the case of Table 2, the solutions of dynamic 
constraints follow from the similarity between force output (fT6l) and elastic-force patterns. 
The resulting dynamic states are discussed in terms of the scaling exponents for the muscle 
dynamic structure [A rm , L rm ] preserving muscle shape and volume ©• Other related ob- 
servable mechanical characteristics are exemplified in Tables 1 and 2. The major outcome 
of the analysis in Tables 2 and 3 is that both slow-twitch and fast-twitch fibres belonging to 
the same muscle m should manifest concerted behavior coordinated by the dynamic active 
elastic forces. 

Another significant feature of the analysis in Table 3 is a striking prediction of the me- 
chanical functions which are expected to be shown by a given striated muscle m of certain 
specialization {primary functions indicated by regime r = 2, see proof below) when its cy- 
cling dynamics is switched to regimes r = 1 and 3 by tuning to the corresponding natural 
frequencies T~^. In case of regime r = 1, both types of arbitrary slow muscle tuned to T± slow 
and fast muscle tuned to T^} ast (Table 4) are expected to show maximum workloop efficiency 
when acting as controlled spring. In the efficient nonlinear regime r = 3 the slow and fast 
struts (m = 3 in Table 2) will not show another function, but any type of brakes (m = 2 
in Table 2) will work as motors, whereas motor are expected to expose a new function, say, 
m = 6 [determined by t]q = (6/7)/(l/7) = 6] that is closer to the brake activity (i]2 = 3) 
than the strut (7/3 = 00). The cardiac muscles seem to display a crucial dynamic state, 
say m = with r] = 0, which flatters the heart. Such predicted secondary functions and 
unusual (m = and 6) muscles adapted to new functions is a challenge deserving further 
study by experimentalists. 

B. Direct observation of muscle specialization 

"If a muscle is specialized for a particular mechanical role how this is reflected in it 
architecture?" [40]. The stated problem is approached here by the comparative analysis 
between the muscle allometric exponents and those predicted for particular efficient activities 
describing the trends of biomass accommodation via PCSA and along a muscle. 
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1. Isolated muscles in hindlimb of mammals and birds 

In Table 5, the morphometric data on the allometric exponents for the mean cross- 
sectional area and length of four skeletal muscles in the mammalian hindlimb 
for 35 quadrupedal species of body-mass domain exceeding four orders in magnitude are 
studied. 

Place Table 5 



First, let us verify the cylinder- shape similarity of skeletal muscles described by Eq. (Q. 
The muscle mass index a 0m estimated in Eq. (fTUj) via experimental data aorn and is 
compared in Table 5 with the measured indexes a^^\ 

Place Fig. 2 



In Figs. 2 and 3, the method of determination of the primary mechanical function is 
illustrated: the adapted muscle structure is indicated by the appropriate theoretical point 
located most closely to the datapoint. 

Place Fig. 3 

The found reliable estimates were used then in the muscle-function analysis in Figs. 
2 and 3. The established small indices aQ m generally validate the muscle biomechanics by 
proving a high-precision observation of locomotory muscle patterns via muscle morphometry 
and functional physiology. This implies that the effect of biomechanical adaptation of muscle 
design to active elastic forces predominates over effects of biological adaptation assigned to 
small «om 

Secondly, the analysis in Figs. 2 and 3 indicates strong correlations between the mor- 
phometrically characterized structure of skeletal muscle and one of the primary locomotor 
functions described in Table 2. The primary functions indicated in Table 5 are found with 
a high degree of certainty. Indeed, as illustrated in Fig. 2, the deviations of distances 
measured along the dashed line, corresponding to a given muscle, between the datapoint 
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and distant challengers for the primary function, from the smallest distance indicating the 
primary candidate, always exceed the experimental uncertainty. 

Thirdly, the found muscle mechanical specifications do not conflict with the physiological 
categorization established for joint extensors and flexors, which muscle structures are shown 
to be adapted to the brake and motor functions via activation of eccentric and concentric 
elastic forces. The found structure parameter r} p i ant ~ 18 indicates the foot support activity 
for plantaris as the primary function (Table 5) that is in accord with in vivo workloop 
presented in Fig. 1C. As shown in Table 3, the struts are most conservative muscles no 
changing their support function in non-linear regimes. In contrast, the gastrocnemius in 
mammals manifests their motor, strut, and brake functions in, respectively, uphill, level, 
and incline running of animals. Through the motor adapted structure with i] gast « r\\ — 4, 
the analysis in Fig. 3 establishes the motor activity for gastrocnemius as the primary 
function naturally selected for the significant mechanical task of uphill running exploring 
the bilinear muscle dynamics. The effective trilinear gastrocnemius-displacement dynamics 
is most close to the brake-like activity rje = 6, attributed to the secondary function of the 
motor experimentally observed in gastrocnemius of incline running turkey [27] and hopping 
tammar wallabies [24]. 

In Fig. 4, the overall muscle peak stress data measured in limb muscles of animals in 
strenuous activity, reviewed by Biewener [25], are re-examined and re-analyzed accounting 
for the primary functions of hindlimb muscles established in Table 5. 

Place Fig. 4 

The uphill-motor specialization of gastrocnemius is independently supported by the 
compressive-stress analysis made in Fig. 4 for fast running, jumping, and hopping mam- 
mals. The stress scaling exponents (s m ) predicted for the motor (si = 1/5), strut (S3 = 0), 
and control (54 = 0) functions are shown to be distinguishable in work-specific mammalian 
muscles described in Table 2. Hence, although the overall-function data by Biewener [25] 
indeed expose almost weight-independent muscle stress, earlier postulated by McMahon in 
Eq. (j3j) and only in part justified here by the slow-fibre muscles (Table 1) and strut muscles 
(Table 2), the analyses in Fig. 4 demonstrates how the function-specific muscle stress may 
serve as a new tool for the direct observation of muscle specialization ignored in all previous 
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overall-function analyses. 

I have also investigated an interesting question: whether the primary function established 
for a certain leg muscle in mammals specialized to fast running coincides with that for the 
same muscle in birds? The pioneering data on individual leg muscles in 8 running birds, 
ranging in size from 0.1 kg quail to 40 kg ostrich, are analyzed in Table 6 and Fig. 5. 

Place Table 6 

Place Fig. 5 

In running and non-running birds (Fig. 5), the gastrocnemius is employed as the brake 
and spring, in contrast to the motor function in mammals (Table 5). This is in accord with 
Bennett [23], who noted that "the full force-generated capacity of gastrocnemius is only 
used occasionally, such as during take-off, when a bird attempts to throw itself into the 
air". This explains our indirect observation: the primary function of the gastrocnemius in 
running specialists is attributed to the foot flexor in mammals and ankle extensor in birds 
(Table 6). In non-running birds, the legs are designed to control the ground locomotion 
(Fig. 5), whereas the wings may share motor and brake functions (Table 3), in accord with 
the review by Dickinson et al. [6]. 

2. Micro-macro scale correspondence 

There are many striking examples when skeletal muscles expose adaptation to a specific 
function, e.g. [43, 3]. The striated muscles anatomically suited to concentric or eccentric 
work [2] are structurally distinct having, respectively, long thin cells or short wide cells [5]. 
This observation suggests the microscopic level of muscle-cell adaptation introduced here by 

A {conc) decent) , j (ecent) j (cone) / 9Q x 
A cell > A cell aI1Cl L cell > L cell \ Z ^) 

for the cellular cross-sectional area A ce u (= A s ) and length L ce u (= L s ). Adopting these 
function specific trends, one may expect to observe the cell-structure parameters n s = 4 and 
3 for sarcomeres accommodated to efficient shortening or stretching of muscle as a whole. 
A general question arises whether allometric coefficients of proportionality omitted above 
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in all structure-function power-law (scaling) relations are also attributed to active elastic 
strains accompanying maximum force production? Or, alternatively, other microscopically 
justified mechanisms, c.f. [44], or additional parameters (such as pinnate angle) may result 
in different general macroscopic consequences? Given the highly conservative nature of 
contractive units of skeletal muscles and their well pronounced organization [25], the specific- 
function trends of the muscle cross-sectional area 

A(isom) a (cone) Aeccen) Asprin) /oq\ 
■^strut -^motor brake -^contr 



and muscle-fibre length 



j (sprin) j (eccen) j (cone) j (isom) (Ql\ 
L contr > '-'brake > L motor > L strut 



are generally expected from Table 2. The suggested trends become observable via the 
primary functions established in Table 5 for gastrocnemius (m = 1), DDF (m = 1), CDE 
(m = 2), and plantaris (m = 3), when the regression data [22] on passive-muscle structure 
[A£ p) (M), L { ^\m)} are taken additionally into consideration: A^ t > > A^ ] F > 

A%*1 and > l£? > « > start ing with M > 1 kg. 

Similarly, the trend for active stiffness 

> > and, generally, > (32) 

straightforwardly follows from Table 2. Given that the optimum velocity for fast fibres 
Vim oc lU 2 (Table 1), Eq. @Q provides 

y(°Vt) y(opt) y(opt) /oo\ 
brake 'motor " strut V' 3 '-'/ 

Moreover, a crude estimate for the cost energy 

CutT2 > CU%$ > CU%£ (34) 

follows from CUf*^ oc M m (122]) and the experimental data by Pollock and Shadwick [22], 
M 1 ^ exp ' ) > Mg 05 ^ > M^ xp \ considered at the same body mass M. The finding (I34|) is 
in accord with the experimental observation [44]: muscles contracting nearly isometrically 
(strut function) generate force more economically than muscles involved in concentric work 
(via motor function). 
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3. Muscle dynamics of mammalian legs and dragonfly wings 



Given that mammalian leg extensors are active mostly during lengthening [2], the brake 
primary function (m = 2 in Table 2) could be assigned to leg muscles specified by effective 
length Li eg oc M 1 / 4 (cti eg = is adopted). In accord with Hill's second constraint, underlaid 
by the proper frequency T 3m ~ oc L m 2 (Table 4), the theory predicts vj^ 3 ^ oc L^ eg oc M -1 / 4 
that results in l/T/™ ax) a L~ 2 g oc M" 1 / 16 . Similarly for the wing-motor muscles in fly- 
ing birds (m — 1 in Table 2) one should expect V^™^ cc L~] ng oc M -1 / 5 , for contraction 
velocity, and oc ikT 1 / 25 , for the frequency or, alternatively, l/T^ ] g oc M" 1 / 5 , in 

the optimum- velocity regime (see Table 4). Hence, analytically revealed Hill's constraint 
becomes observable via the empirical regression data by Medler [43]: on the maximum- 
amplitude contraction velocities for the locomotor muscles in leg of terrestrial animals, 
Vieg^ 0C M~ ' 25 , and that for wings in flying birds, bats, and insects, V^^J oc M~ a2 °. 
Moreover, the experimental data by Schilder and Marden [45] of the wingbeat frequency 
1 lowing o^ M m - 20 scaled by mass M m (and length L 0m ) of the basalar muscle in dragonflies 
indicate that the motor-type muscles (L 0m cc Mm 5 , see analysis in Fig. 6) were studied 
self-consistently in the optimum, steady-velocity motion regime. 

Place Fig. 6 

In the same optimum- velocity regime (Table 1), the maximum-amplitude static force 
AF^ ow) oc M™ /3 and net lever-system force F^J ] AF^ st) oc M m reported 
by Schilder and Marden [45] may be associated with the slow and fast activated fibres 
in the basalar muscles tuned elastically to the linear regime through the dynamic PCSA 
^tnf 1 ^ Mm 3 and length cc mU 3 ■ The observed dynamic force output Ff g ^ oc M m m 

can be therefore suggested as the mixed-fibre force F^J^ = AF^™ 1 ^ cc Mm 6 (Table 1), i.e. 
as a compromise of the forces F^t^ an d F^ p \ These estimates challenge further analysis 
of the reported dynamic forces. 
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V. CONCLUSION 



A theoretical framework for mechanical characterization of the three transient activated 
states of the striated muscles passing in force-length cycles through the three distinct dy- 
namic regimes is proposed. The explicit analytical description of muscle locomotor func- 
tions and related mechanical characteristics is provided on the basis of two concepts: (i) the 
preservation of spindle-type shape in skeletal muscles and egg-type shape in cardiac muscles 
related to the preservation of dynamic muscle volume and (ii) the mechanical similarity be- 
tween action and reaction forces emerging in biomechanically equivalent states. Exploring 
known patterns of elastic forces in continuum mechanics, the macroscopic study of the force 
production and its functional and structural accommodation in the loaded muscle organ as 
a whole provides the following major points. 

1. It is demonstrated how the dynamic (frequency- velocity) constraints for muscle con- 
tractions, first observed by Hill in hovering birds and then revealed in locomotor muscles of 
running animals and flying birds, bats, and insects, can be derived from the generic principle 
of mechanical (force and velocity) similarity. 

2. It is shown how relations in classical mechanics of solids can be explored in soft 
tissues. The study is grounded by the active-force muscle stiffness reliably derived in all 
muscle work loops nearby and below the maximum-amplitude exerted forces. The muscle 
stiffness, underlaid by sarcomere stiffness, is shown to be dependent on muscle geometry and 
dynamic functional variable, underlaid by elastic moduli, which encompass all contractive 
elements acting as an elastic continuum medium. 

3. The theoretical prediction that the fast and slow muscles should generate maximum 
forces linear, respectively, with the muscle volume and cross-sectional area, regardless of 
muscle function and structure, is in part validated by the direct empirical observation of 
maximal forces exerted by animals and by the provided indirect observation of the adapted 
(primary) muscle functions in legs of mammals and birds. 

4. The macroscopic structures of locomotor skeletal muscles observable directly by mus- 
cle allometry are found to be adapted to the maximum-force state, following moderate- 
velocity dynamic regime, instead of the expected optimum velocity regime. Such a bilinear- 
displacement muscle dynamics involving both fast-twitch and slow-twitch powering muscle 
fibres sheds light on the origin of allometric power laws and muscle specialization. The 
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adapted structures are examined via available empirical data: the legs are brakes in mam- 
mals and springs in non-running birds, whereas the wings are motor-brake engines in flying 
species. Suggested pump function for the cardiac muscles needs further experimental tests. 

5. The provided study of the muscle specialization in mammalian hindlimb indicates that 
the properly tuned force production function is a dominated factor in the accommodation 
of muscle structure. This finding also indicates the predomination role of mechanical effects 
over biological adaptive mechanisms assigned to the relatively small muscle-mass index. As 
the result, a new investigation tool for indirect statistical observation of the biomechanical 
adaptation of individual locomotor muscles is proposed through the regression analysis of 
in vivo muscle stresses in synergists scaled across different-sized animals. 

6. The assumption on that the muscle tuning muscle ability of animals can be modeled 
by active elastic forces via non-linear muscle elastic moduli is validated by the observation 
of the theoretical predictions for muscle dynamics of legs and wings in running and flying 
specialists. Predictions are made for the experimental modelling the primary and secondary 
function by tuning the cycling muscle to the corresponding natural frequency and controlling 
its efficiency. 

7. The conservative character of architecture and related mechanical characteristics of 
striated muscles suggests general trends following from mechanical and shape constraints. 
The trends dictated by primary functions explain, in particular, why the muscles having 
larger fibre and sarcomere lengths and suited to efficient eccentric work, tend toward higher 
optimum contraction velocities, but show lower maximum stiffness and mechanical energy 
cost. 

8. As an intriguing outcome of the analysis of maximal contraction muscle velocities 
and frequencies, the maximum-speed steady locomotion is revealed to be controlled by non- 
linear elasticity of slow-fibre muscles generating moderated force. This finding deserves 
further evaluation in finite muscle element analysis studying top speeds of living and extant 
animals. 
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Appendix A. List of abbreviations 

PCSA - physiologic cross-sectional area 

Mathematical signs and symbols 

= - common equality sign 

= - identity sign implying "by definition" 

~ - approximate equality sign 

~ - proportionality relation symbol omitting only numerical coefficients 
= - here used as similarity sign supporting only physical dimension units 
oc - here used as scaling rule symbol not supporting dimension units 

Physical and geometrical notations 

a m - muscle-mass allometric index 

em Pt ^- muscle strain in the optimum dynamic regime 

i] m - muscle geometry parameter 

puss- tissue density 

a uss^~ P e& k tissue stress 

AL - length change 

AF - force change 

At m - activation timing of muscle m 

A rm - cross-sectional area of muscle m in passive (r = 0) and active (r ^ 0) states 
a- scaling exponent for cross-sectional area 
D - diameter of ideal cylinder 

E rm - active-muscle elastic modulus establishing the dynamic regime r ^ 
e- strain scaling exponent 

AF^f* = AFm* 1 ^- maximum active elastic force 
Fprod^' production force by fast muscle 

Fmotor' elastic force adapted to concentric work in motor muscle 
i 7 m«sc' ) - maximum force exerted by muscle 
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K m - active muscle stiffness 
K s - sarco mere/cellular stiffness 
L- length of an ideal cylinder 

L m - variable muscle length in non-specified dynamics 

Lrm - dynamic muscle length in the regime r 

/ - length exponent 

m - muscle in unspecified function 

M - body mass of animals 

M m - muscle mass 

r - numerical parameter indicating transient dynamic states via optimum- velocity (r = 1), 
moderate- velocity (r = 2), and high- velocity (r = 3) dynamic regimes, distinct of passive 
muscle state (r = 0). 

T rm - period of cycling in the adapted regime r 

V rm - muscle contraction velocity in the dynamic regime r 

W - body weight 

Appendix B. Scaling Muscle Functions 

The motor function is associated with the active force generated during muscle 

shortening at moderate contraction velocity at the turning points 2 in Figs. 1A, IB, and 
1C. In fast-fibre muscles, the corresponding concentric force 

rp(conc) _ r?(m&x) _ » pfomc) ci(/ast) , 3/2 T -\ jp(fast) Mr^ 
r motor — r elast ~ LAr 2m Cj 2m I± 2m- U 2m — r prod \°° J 

is described by the known universal pattern of the maximal elastic forces [33] equally applied 
to pure bending, pure torsion, or complex bending-torsion loads subjected to long cylinder 
of length L 2m and cross-sectional area [26] A 2m - The exploration of Eq. f)35p though Eqs. 
® , COID ) (jUD ) and fl26j) results in the fast-muscle-force constraint 3a m /2 — l m = 1 + a m . It 
is remarkable that the case of slow- fibre muscle, namely 

^(conc) _ r (max) _ A ^(conc) rp(slow) ,3/2 T -1 ^ ^(slow) / q „n 
* motor — * elast ~ ^ 2m ~ &2m A 2m L 2m ~ b prod K 6 ®) 

results in the slow-muscle-force constraint 3a m /2 — 2l m = a m , which is exactly the same as 
fast muscle, in view of function-independent Eq. ffTUl) . Therefore, any muscle tuned to the 
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motor locomotor function should expose its dynamic structure scaled by 

(cone) 4 -.(cone) \(\ \ \ ( r K7\ 

a motor gl- 1 - ' & motor ) > ''motor ' ^motorji I"' J 

regardless of the fibre type content. This finding follows from both the muscle force con- 
straints solved with the help of the function-independent muscle-shape constraint (11 01) . 
Moreover, as shown in [26], the principal component of the compressive stress <Jm° n ^ speci- 
fying Eq. (j!7p may be caused by the peak transverse-tensile strains 

A 7-)( max ) n 
(cone) _ £m _ (cone) _ «m _ , /on\ 

t motor —j ^ lvl m ) wil > Li ^m — ^motor ~ l mi \°°) 

where AD^ ~ .D m ^ Am 2 is transverse muscle deformation. 
Likewise, the maximum elastic eccentric force 

^(eccen) A ^(eccen) ^(eccen) A 2 T -2 ~ ^(fast) ( oc\\ 

? brake - lXt 2m ~ -^2 m ^m-^m = b prod > l dy J 

associated with the 6raA;e muscle function (Fig. IB) provides the maximum elastic stress 

j-i(eccen) 

(max) _ f brake M s m __ (eccen) _ „, /^s 

a 6rafce ~ 1 m ' m ~~ &™fce ~ a m ~ ^m, ( 4 UJ 

following from Eqs. (117j) and ( 1391) . The unique solution to both fast-muscle-force constraint, 
2a m — 21 m = 1 + a m , and slow-muscle-force constraint, 2a m — 3l m = a m , is 

(eccen) 3, , i(eccera) (eccen) /i i \ / j -. \ 

a bmfce - ^l 1 + "6rafceJ , We ~ S ferafce ~ 4 1 1 + a brake ) ■ (41 J 

The strut muscle function treated as antagonistic to both motor and brake functions 
drives nearly isometric contractions characteristic of small, but non-zero length change 
(AL m <C L m ) achieved near peak forces (see Fig. 1C). This suggests the nearly isomet- 
ric force 

T-,(isom) _ rp(isom) (isom) , ^ r-i(fast) -.i (isom) _ A j (isom) / j (A0 \ 
P strut — tj 2m £ 2m A 2m — p prod ) Wlttl £ 2m ~ lXLj 2m / L 2m, 

in fast muscles. Again, one solves the muscle strut constraints 2a m + l m = 1 + a m and 
2a m + 2l m = a m resulting in 

(isom) _ , j j(isom) _ (isom) _ n . , a j (isom) j2 

a strut — 1 i "strut ana ^ trnt — S sirui — U, Willi /A^ 2m CC L 2m , 

for any type of muscles. 
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A new antagonist (to strut muscle) tuned to the cardiac type contractions via active 
elastic force 

rn(card) _ A tp {card) rp(card) T 2 ~ rp(fast) (..s 

fpump ~ A ^2m ~ ^2m L 2m = k prod I 44 J 

is associated with, say, pump function providing the fast-muscle-force constraint 2l m = 
1 + a m . This yields 

(card) _ j(card) _ (card) _ (card) _ \(-> , \ (Ar\ 

U'pump h pump l3 pump '-'pump X 1 - ~ XJL pump) i \^' J J 

equally applied to slow-fibre muscles resulting in the slow-force constraint l m — a m . 

To complete the intrinsic-force description, the spring-type control function associated 
with the optimum-regime elastic force 

= F^l a E[T ] MT « E^A^L^ = (46) 

in slow-fiber muscles results in 

As^n) _ 2 (spnn) _ 1 ( , (sprin) _ ^sprin) _ , . 

a cont — gl 1 + a cont), <- C ont ~ gV 1 + a cont), WltU S cont — e cont — U, (_4 < ) 

that follows from the slow-force and fast-force constraints 2(a m + l m )/3 = a m and (2a m + 
5/ m )/3 = 1 + a m and therefore is valid for any type of muscle tuned to velocity-optimum 
regime. All obtained specific-function mechanical characteristics are summarized in Table 
2. 
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Figure Captions 



Fig. 1. The qualitative analysis of the in vivo muscle force-length data. The muscle 
motor function is presented by gastrocnemius powering during shortening in uphill running 
turkey (inset A, adapted from [27]). The lateral gastrocnemius and plantaris act as brake 
(inset B) and strut (inset C) in hopping tammar wallabies [24]. The solid (and dashed) 
arrows indicate rasing (and decreasing) of the exerted force near its maximum magnitude 
-Pmax- The regions of the linear force- length domain are displayed by the force change AF lm 
and length change AL im , both estimated from L im , and the starting datapoint F\ m of the 
force enhancement achieved at the optimum contraction velocity V opt . Similar to physical 
pendulum, the resting length L 0m is expected to be passed at near maximum velocities V m3x 
and lower forces F^ m . The origin of intrinsic muscle forces (inset D): in both cases of the 
powering shortening (motor) and lengthening (brake) muscles the resulted force -F max is a 
superposition of the production force output F pro d and reaction passive F pass and active F act 
elastic forces [see also text below Eq. (TL6l)]. 

Fig. 2. The indirect observation of the primary activity of mammalian plantaris. The 
solid symbol is the datapoint [22] presented in Table 5 and the bars indicate experimental 
error. The open symbols are theoretical estimates for stable dynamic structures established 
for the motor, brake, strut, or control functions described in Table 2, with a m = taken 
from Table 5. 

Fig. 3. The observation of the primary mechanical function in some isolated indi- 
vidual muscles in mammals. The analysis and notations correspond to those in Fig. 2. 
The experimental (and theoretical) data for gastrocnemius, DDF (deep digital flexor), and 
CDE (common digital extensor) are shown, respectively, by the closed (and open) inverted 
triangles, regular triangles, and circles. All the data are taken from Table 5. 

Fig. 4. The qualitative study of the in vivo data on the peak stress in individual leg 
muscles of animals in strenuous activity. The symbols employed above in Figs. 2 and 3 
are extended by the open circles (triceps) for the data on peak muscle stress taken from 
Table 1 in [25], with the exclusion of the slow-mode data on cantering goat and trotting cat. 
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The data [41] on the activated isometric stress in isolated white rabbit tibialis are added. 
The dashed line shows the brake-functional stress indicated by the stress scaling exponent 
s = 1/4. The solid lines are drawn by 115 • M 1//5 , for the motor function, and by 215 kPa, 
for the strut and spring functions. All coefficients are adjusted by eye. 

Fig. 5. The analysis of the primary mechanical functions for leg muscles in running 
and non-running birds. The measured (and estimated) data taken from Table 6 (and Table 
2) for gastrocnemius, femorotibialis, and digital flexors are shown by the closed (and open) 
inverted triangles, circles, and regular triangles, respectively. The semi-open triangles are 
the data by Bennett [23] for non-running birds. 

Fig. 6. The qualitative scaling of the basalar structure to muscle mass in male dragonflies 
(Odonata and Anisoptera, listed in Fig. 5 in [45]. The datapoints for muscle length is a 
courtesy by the authors. The estimated muscle cross-sectional area is obtained on the 
basis of Eq. (CQ) taken with pmSc = 1060 kg/m 3 [46]. The solid lines are L motor = 0.052 -Mm 5 
and A mo tor = 0.018 • M m . The dashed lines indicated by the scaling exponents are drawn 
according to muscle specialization shown in Table 2. All pre-exponential coefficients are 
adjusted by eye. 
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Tables 



Optimum muscle characteristics, (equations) 


Fast fibres 


Slow fibres 


Mixed fibres 


Optimum length change, ALi m , fl25|) 


Lm 




Lm 


Production/active-elastic force, AFi m , (fl"6j). (Tl8j). ([25]) 


A T 


A 


A T 1 ' 2 

^m lj m 


Optimum stiffness, = Ex m A Xm / L Xm , ([19]) 


A m 


A J- 1 


A T~ 1/2 


Optimum elastic stress, o\ m = AFi m /Ai m , (TTT1) 


L m 


^m 


7-V2 
J-'m 


Contraction frequency, ~ \/E lm / p 0m /L lm , (USD 


T -1/2 
J-'m 


r-1 

^m 


r-3/4 
J-'m 


Optimum velocity, Vj. m = k!Ts*c, (EH]) 


7-1/2 
J-'m 


^m 


7-1/4 
J-'m 


Optimum power, P lm = Fi m V lm 


A r 3/2 


A 


A r 3/ * 

rt-mJ-'m 



Table 1. General mechanical characteristics of the striated muscles tuned to linear- 
displacement dynamic regime scaled to dynamic fiber length L m = L lm . The mixed-fibre 
scaling dynamic exponents (shown in the last column) are modeled by the common means for 
the fast-muscle and slow-muscle exponents (established in the second and third columns), i.e. 



F mix ^ >y Ff ast F s i ow ; A m and L m are attributed to the stabilized dynamic muscle geometry 
constrained by muscle volume ([9]). 
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Locomotor pattern, regime 


Motor, r=2 
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Table 2. The locomotor functions and their mechanical characteristics scaled to dynamic 
structures. The all-type powering individual muscles m = 1,2, 3, and 5 are tuned to the 
maximum- force bilinear dynamic regime r = 2 [described in Eq. (|26p ] and muscles m = 4 
act in the linear regime r = 1 [Eq. ( l25l) ]. *The data shown for fast muscles. The allometric 
exponents are related to animal's body mass via Eq. (jSJ). 
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Table 3. Locomotor functions predicted by dynamic structured for slow and fast striated 
muscles tuned to distinct dynamic regimes. The primary functions (r = 2) are shown by bold 
type. The analysis of functional muscle structures made in terms of elastic-force patterns: 
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the active- muscle optimum- velocity (r = 1), moderate- velocity (r = 2), and high- velocity 
(r = 3) dynamic regimes are described in the first column via the muscle elastic moduli E rm 
[Eqs. ([25]), (126")) . and (127)) ] and specified by slow and fast force output [Eq. (J2JJ], shown in 
the second column. The third and next odd columns show the elastic force functional scaling 
in concentric, eccentric, isometric, and pump contractions. The corresponding solutions to 
scaling equations underlaid by the force similarity principle ( 1161) are shown for simplicity 
with a rm = 0, in the forth and next even columns. Notation: nc indicates non-conclusive 
solution. 



Dynamic regimes 


Optimum, r = 1 


Moderate, r = 2 


Maximum, r = 3 


Muscle type 


slow 


fast 


slow 


fast 


slow 


fast 


Natural frequency, Eq. (fl5l) 


r-l 

m 


r - 1/2 


r-3/2 


r-i 


J- 2 


r-3/2 
J-'m 


Contraction velocity, Eq. (1231) 


r o 
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r- 1/2 


ro 
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Table 4. Dynamic characterization of the red (slow) and white (fast) striated muscles in 
the optimum-, moderate-, and maximum- velocity dynamic regimes r = 1, 2, and 3 described 
in Table 3. 



Individual mammalian muscles 


(exp) 


,(exp) 
L 0m 


(exp) 


V0m 


(est) 




4n 


Prim, functions 


Gastrocnemius (and soleus) 


0.77±.02 


0.21 ± .02 


-0.03 


3.7 


-0.02 


0.78 


0.20 


motor, m = 1 


Deep digital flexor (DDF)* } 


0.85 ± .03 


0.18 ± .02 


0.03 


4.7 


0.03 


0.82 


0.21 


motor, m = 1 


Comm. digit, extensor (CDE) 


0.69 ± .04 


0.24 ± .02 


-0.07 


2.9 


-0.07 


0.70 


0.23 


brake, m = 2 


Plantaris (SDF) 


0.91 ± .04 


0.05 ± .04 


-0.03 


18 


-0.04 


0.96 


0.00 


strut, m = 3 


Ankle-joint muscle group 


0.81 ± .03 


0.17 ± .03 


-0.03 


4.8 


-0.03 


0.78 


0.19 


motor, g = 1 



Table 5. The analysis of the allometric data by Pollock and Shadwick [22] provided 
on the basis of Eq. (ITOj) and Table 2. The shown statistical error is approximated by the 
symmetrized 95% confidence interval. The methodology of the analysis is illustrated in Fig. 
2. The primary functions found in Figs. 2 and 3 are described following Table 2, with 
. The overall muscle group (g = 1) is determined as the standard mean over all 



a.. 



a 



(est) 
Om 



muscles. *^DDF includes individual flexor hallucis and flexor digitorum longus; SDF means 
superficial digital flexor. 
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Running birds 


„( cx p) 


(exp) 


j(est) 


„(exp) //(est) 
0m 1 0m 


0>2m 


hm 


Primary function/force 


Gastrocnemius 


0.81 ±0.14 


0.14 


0.33 


2.5 


0.85 


0.29 


brake / eccentric 


Digital flexors (DF) 


0.76 ± 0.22 


-0.03 


0.21 


3.6 


0.78 


0.19 


motor /concentric 


Femorotibialis 


0.80 ±0.12 


-0.02 


0.18 


4.4 


0.78 


0.20 


motor /concentric 


Overall group 


0.79 ±0.16 


0.03 


0.24 


3.3 


0.77 


0.26 


brake/eccentric 



Table 6. The analysis of the allometric data by Maloiy et al. [42]. The shown large 
error is due to relatively wide confidence limits. The mean exponents l^f 1 are estimated via 
Eq. ([TO]) . The overall muscle group is determined as the standard mean over all muscles. 
The indicated primary functions and active elastic forces are described by the evaluated 
dynamic-structure exponents aim and li m found as most close to the experimental resting- 
volume data on a e ^ and l£ p ' and therefore assigned to regime r = 2 (Table 2). 
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This figure "Slidel.GIF" is available in "GIF" format from: 



http://arxiv.org/ps/0909. 1444vl 



This figure "Slide2.GIF" is available in "GIF" format from: 



http://arxiv.org/ps/0909. 1444vl 
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